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ABSTRACT 


This paper describes a numerical study of the g-jitter driven double diffusive convective 
flows, thermal and concentration distributions in binary alloy melt systems subject to an external 
magnetic field. The study is based on the finite element solution of transient 
magnetohydrodynamic equations governing the momentum, thermal and solutal transport in the 
melt pool. Numerical simulations are conducted using the synthesized single- and multi- 
frequency g-jitter as well as the real g-jitter data taken during space flights with or without an 
applied magnetic field. It is found that for the conditions studied, the main melt flow follows 
approximately a linear superposition of velocity components induced by individual g-jitter 
components, regardless of whether a magnetic field exists or not. The flow field is characterized 
by a recirculating double diffusive convection loop oscillating in time with a defined frequency 
equal to that of the driving g-jitter force. An applied magnetic field has little effect on the 
oscillating recirculating pattern, except around the moment in time when the flow reverses its 
direction. The field has no effect on the oscillation period, but it changes the phase angle. It is 
very effective in suppressing the flow intensity and produces a notable reduction of the solutal 
striation and time fluctuations in the melt. For a given magnetic field strength, the magnetic 
damping effect is more pronounced on the velocity associated with the largest g-jitter component 
present and/or the g-jitter spiking peaks. A stronger magnetic field is more effective in 
suppressing the melt convection and also is more helpful in bringing the convection in phase 
with the g-jitter driving force. The applied field is particularly useful in suppressing the effect of 
real g-jitter spikes on both flow and solutal distributions. With appropriately selected magnetic 
fields, the convective flows caused by g-jitter can be reduced sufficiently and diffusion dominant, 
solutal transport in the melt is possible. 
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applied magnetic field, its magnitude 
heat capacity 

maximum and minimum concentration 
concentration difference 
dimensionless concentration 
unit vector 

dimensionless frequency of the acceleration 
frequency of the acceleration 
dimensionless amplitude of the acceleration 
amplitude of the acceleration 
earth gravity 

Hartmann number, Ha = BL-Jaj ~Ji 

solutal Grashof number, Gr$ = gficACL 3 /V 2 

thermal Grashof number, Gr T = gfij-ATV /V 2 

unit vectors of the ith, jth components 
thermal conductivity 
length of the square cavity 
discretized nodal pressure array 

* T V 

Prandtl number. Pr = — 

a 

y 

schmidt number, Sc = — 

D 

dimensionless temperature 
maximum and minimum temperatures 

temperature between T H and T L 
dimensionless velocity 

x and y component of the dimensionless velocity 
volume of the cavity 
dimensionless coordinates 


Greek 

a thermal diffusivity 

P c solutal expansion coefficient 

l 3 t thermal expansion coefficient 

D solutal diffusivity 

<f> shape function for velocity 

0 shape function for temperature and concentration 

y/ shape function lor pressusre 



so ^ 


M 

v 

dQ 
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dynamic viscosity 

kinematic viscosity 

boundary of computational domain 

gradient operator 

density 

computational domain 


Subscripts 


i,j the ithjth point 

n the nth component 

Superscripts 

i,j the ithjth component 

T matrix transpose 
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1. INTRODUCTION 



During space flights, perturbation from idealized .dynamic weightlessness conditions 
established in Space Shuttles arises from a variety of disturbances such as crew motions, 
mechanical vibrations (pumps, motors, excitations of natural frequencies of spacecraft 
structures), spacecraft maneuvers and attitude, atmosphere drag and the Earth's gravity gradient 
[ 13 - Some typical signatures of gravity perturbation or more often referred to as g-jitter, taken by 
astronauts using on-board accelerometer during a typical Space Shuttle flight, are illustrated in 
Figure 1^ Although small in magnitude, g-jitter is random in both direction and time, and causes 
time dependent convection in melt being processed in space. Experiments on melt growth of 
crystals conducted in Space Shuttles suggest that the g-jitter induced melt flow has deleterious 
effects on the solute uniformity and defect^ formation in the space grown crystals and is the main 
cause for the non-repeatability of the crystal quality processed in microgravity [2,3], 


A 


Direct study of the g-jitter effects on the thermal and fluid flow behavior of melt growth 
systems in space environment is difficult because of limited instrumentation precision. 
Numerical simulations paw^^-be: useful in providing a basic understanding of the physics 
governing the g-jitter induced phenomena. Consequently, both 2D and 3D numerical models 
have been developed and applied to estimate the effects associated with both idealized single- 
and multiple-frequency g-jitter modulations and realistic g-jitter data collected by accelerometer 
during actual flight experiments. Calculations suggested that the frequency, amplitude and 
spatial orientation of the residual gravity vector all play an important role in determining the 
convective flow behavior of the system. The most detrimental effects occur when g-jitter is 
perpendicular to the applied thermal gradient [4-7], For solidification or crystal growth 
applications, space experiments and numerical simulations both indicate that g-jitter components 
with lower frequencies are most harmful and are repsponsible for defects and solutal striation in 
the crystals produced in Space Shuttles. 


The need to have a quiescent melt pool from which crystals are grown has motivated 
researchers to develop and apply various concepts to suppress the deleterious g-jitter induced 
convection. Because of its proven success in reducing thermally induced convection during melt 
growth of semiconductor crystals under terrestrial conditions, magnetic damping has been 
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considered as a potential technique to eliminate the harmful g-jitter convection, or more correctly 
stated, to damp the convection beyond detection. Some basic studies have been carried out in 
this direction and numerical simulations have been conducted also for the melt systems where 
flow is primarily attributed mainly to thermal gradients [8,9], 

Thus far, very scarce work appears to have been done on the fundamental understanding of 
magnetic field effects on double diffusive convection in microgravity environment. The study of 
this type is of critical importance to space experiments involving binary alloy solidification. 
Some of these experiments are being planned for future Space Shuttle flights [7]. Similar 
experiments are being considered for the upcoming International Space Station. In this paper, 
we intend to develop a fundamental understanding of magnetic field effects on double diffusive 
convection induced by g-jitter in microgravity environment through extensive numerical 
simulations and analysis. For this purpose, a simple 2-D cavity geometry is considered and the 
Navier-Stokes equations, energy equations and solutal transport equations formulated within the 
framework of magnetohydrodynamics are solved using the finite element method. Numerical 
simulations are performed with the g-jitter data synthesized using Fourier series transformations 
and real g-jitter data taken directly from space flights, in combination with an external magnetic 
field. Computed results suggest that an external magnetic field can have a significant effect on 
double diffusive convection and a quiescent melt pool is possible with a strong field. The 
information and computational methodology should be useful for an understanding of physics 
governing the g-jitter induced double diffusive convection in a magnetic field and for designing 
and operating thermal fluids systems for space experiments involving electrically conducting 
binary melts. 


2. PROBLEM FORMULATION 


The model problem under consideration is schematically illustrated in Figure 2, where the 
coordinate system for analysis is also shown. A square cavity is filled with a binary liquid metal 
alloy melt. Its right side is fixed at a higher concentration C H and its left side a higher 
temperature T H . One of the important characteristics is that the solutal Grashof number Gr s and 
the thermal Grashof number Gr-r to be approximately the same. The dimension of the square is 
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set to L, as nondimensionalized, equal to 1 . The upper and lower walls are kept adiabatic. The 
origin of the coordinate system is fixed at the lower left comer of the cavity. In microgravity, 
double diffusive melt convection arises because of the combined effect of the temperature and 
concentration gradients and the time varying gravity perturbation. The melt flow in turn will 
alter, in principle, both the thermal and concentration fields. A DC magnetic field is applied in 
the x-direction. It interacts with the velocity field to produce an opposing Lorentz force to the 
local velocity, thereby resulting in the local velocity being reduced. Consequently, the 
convective transport of solutal elements is supressed. 

In general, g-jitter is known to be random both in time and in direction (see Figure 1). 
However, for a certain period of time, g-jitter may be partially represented by a synthesized 
Fourier series via the fast Fourier transformation, 

g*(f*)= Tgn sin(27z/,*/*)c (1) 

fl — \ 

where g n ’ and f n * are the amplitude and frequency of the nth acceleration component, 
respectively. Here the subscript * designates a dimensional quality. 


The magntohydrodynamic form of the transport equations describing the g-jitter induced 
melt flow, heat and mass transport in the liquid pool system with an imposed magnetic field has- 
been established. These equations are written in a nondimensionalized forms as follows [10]: 


V-u = 0 (2) 

— + (u ■ V)u = -VP + V 2 u - (i Gr r T + Gr s C) g(t) + Ha 2 (u x b x b) (3) 


ST 

dt 


+ U .V7’ = — V 2 T 
Pr 


— + u*VC = — V 2 C 
dt Sc 


(4) 

( 5 ) 
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In the above equations, lengths are nondimensionalized by L, velocity by v/L, time by L 2 lv , 
temperature by AT = T H - T, , and concentration by AC = C H -C L . Also use has been made of 
the Boussinesq’s approximation, p(T,C) = /? 0 (l- /3-p (T- Tq)~ Pc(C-Cq)), and g(t) is 
nondimensionalized by go=9.8 m/s 2 . 

The last term on the right hand side of Eq. (3) represents the Loretnz forces that result 
from the interaction of the applied magnetic field^and the local velocity. These forces always 
act in the opposite direction to the velocity and thus generate a local impedance to the flow. In 
writing Eq. (3) some simplifications regarding the behavior of the magnetic and electric fields 
have been incorporated. First, the effect of the magnetic field induced by the velocity is 
neglected, in comparison with the imposed magnetic field. This is permissible in that the 
magnetic Reynolds number is sufficiently small for the system being studied [11], The theory of 
magnetohydrodynamics has further established that this condition holds true for almost all the 
laboratory experiments involving conducting fluids [10]. Also, with the magnetic field applied 
in the x- and/or y-directions, the electric field (thus the electric potential) arising from the 
interaction of the magnetic field and the velocity can be shown to be zero and thus make no 
contribution to the Lorentz forces [12], 

The above equations are subject to the physical constraints that are described by the 
following boundary conditions: 


u=0 at x=0 & x=l and y=0 & y=l; 

5T dC n p , 

— = — = 0 at v=0 & y= 1 : 

dy dy 

T=l, C=0 at x=0: 

T=0, C=1 at x=l . 

3. THE FINITE ELEMENT SOLUTION 

The above governing equations along with the boundary conditions are solved using the 
Galerkin finite element method. Following the procedures given in [13], the integral forms of 
the governing equations can be obtained with the result, 
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( 6 ) 


(.\nW'-V<t> T dV)V j = -s(\ n y/^ T dV)Y 

( £ <f>fdV) ^. + (^u- Vr)u ; - ( J7 • v<*y r dvyp 
+ (J^V^ • V* : r ^)U / + (J Q (/ • V (Z S)(; • Vf T )dV)Vj + (J n G>*S/*0 7 '</F)T 

+ (Jq Gr$gj<p 0 T dV) C = J^q n • r • / + Jq Her (u x b x b)fidV 


(j^00 T dV)^- + (^0 u • V0 r dV)T + (^~-V0- V0 r dV)T =- J^&iS 


( 7 ) 

( 8 ) 


(^ r ^)^ + (J i ^u-7^ 7 'i/F)C + (| i ^V^V^ 7 '^)C = - i [ Q g e ^ (9) 

The integration defined in the above equations can be carried out numerically over each of 
discretized finite elements once the shape functions <|>, 0, and vj / are appropriately selected. The 
final results may be summarized in the following global matrix, 


"M 

0 

0 ’ 

‘u" 


A(U) + K + jEM^E 7 " B r B c 

"u" 


' F ' 

0 

N r 

0 

f 

+ 

0 D r (U) + L r 0 

T 

= 

G-p 

0 

0 

N C _ 

c 


0 0 D C (U) + L C 

C 

1 

_ G c_ 


(14) 


where the coefficient matrices are computed by 


M P =\ a VV T dV, 

M = f ’ a 4t T dV-. 

Ej-p-V^'dV 

L, = l — VOVO’dd- 

T hp r 


N r = j Q 00 r dV 
N r = [ i OO T dV 

L. = f— V0-V0 T dV 
c *>Sc 

A(U)= -V0 T dV 
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D 7 .(U)= J^u -Vd T dV 
G C =_ I5Q <lc edS 


d c (U) = ^eu-ve T dV\ 

B r = \^Gr T g<f>d T dV 

G T = ~ldQ c jt & cIS ; 

F/ = J^n-r-Z^S 

Ky = (J Q Vf + \n(i-^)U-^ T )dV-{\ a b i bjW T )dV{\-Z >ij ) 

+ {\^b 2 2 b v M T dV)b ij +{\ Q b l 2 b 2 j ^ T dV)b ij 

It is noted that in developing the above matrix equation, the pressure is approximated by the 
penalty method, and P in the momentum equation is substituted byjM^E 7 ^. The assembled 

global matrix equations are stored in the skyline form and solved using the Gaussian elimination 
method. The time integration employs the backward Euler (or implicit) time difference scheme, 
with automatic time step control [13]. 



4. RESULTS AND DISCUSSION 

Extensive numerical simulations were performed using a finite element code that 
incorporates the transient finite element formulation described above. The code was developed 
by our research group over the past years and has been tested against many bench-mark 
problems and other available commercial thermal fluid flow computational packages. It has also 
been used to study the g- jitter induced double diffusive convective flows in the absence of a 
magnetic field [14]. As a result, only the magnetic field effect needs to be added to the code. To 
ensure its predictability of magnetic damping phenomena, the code is also tested against the 
analytical solution obtained for the thermally induced flows in a simple parallel plate 
configuration with a transverse magnetic field imposed [15]. There exists gratifying agreement 


between the numerical and analytical solutions as evidenced by the results appearing in Figure 3, ^ 

thereby providing a validation of the finite element computational methodology presented, v 


It is well known that a Hartmann boundary layer exists near the solid walls whenever an 


external magnetic field is imposed on a moving fluid, and for an oscillating flow the layer 
becomes even thinner [11]. This tiling boundary layer can cause considerable numerical 
difficulties and introduce numerical errors if mesh is not carefully designed. Numerical tests 
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were carried out to determine the grid dependence and computational accuracy «&Ehthe testing 
used the worst scenario condition to be encountered in the study. The testing results suggest that 
about 3 to 5 nodes within the Hartmann boundary layer are required to ensure an accurate 
representation of the wall layer. Combined with the testing of mesh arrangement for 
computation of pure double diffusive convection [14], the final finite element mesh was chosen 
for the computation presented in this study. The mesh consists of 225 9-node elements with 
node spacing decreased from center towards the walls following an arithmetic progressive ratio 


of 0.3. ~ 

For the results presented below, the parameters used for computations are: Pr=0.012,(Gr ~ 0 

^^L54jtH0j_Gr = - 1.54 jcjftj Sc = 65.6 and Ha = 40 ~ 200 unless otherwise indicated. These 
parameters were obtained based on a system being considered for a Space Shuttle flight [14], 

The initial condition for the transient simulations is the steadyrstate fluid flow and concentration 
distribution in the same geometry, without an external magnetic field, which is induced by a 
microgravity field of g=l x 10' 6 go, oriented in the y-direction. This would correspond to a 
typical quiet condition during a space flight. The flow, thermal and solutal distributions for the 
steady state condition., were presented in the previous publication [14] and thus are omitted here. 

It is worth noting, though, that for the low Prandtl number fluids being considered here the 
thermal fields remain basically the same when the g-jitter perturbations become active. For this 
very reason, the thermal contours are not included in discussion. Also, the results discussed 
below are computed based on the situations in which the magnetic field is turned on when g-jitter 
sets in to investigate the magnetic damping effects. This is the optimal scenario from the 
standpoint of power consumption in space. Additional numerical simulations were also 
performed to investigate the cases where the magnetic field is on all the time including the quiet 
period and the results show that the magnetic damping effect is even more pronounced. 

Nonetheless, the results discussed below should serve well the objective of this paper, that is, to 
develop a fundamental understanding of physics governing the effects of magnetic field on the g- 

A 

jitter induced double diffusive convection. 


A. Single frequency g-jitter 


As stated earlier, g-jitter is random in nature and thus a single frequency representation is 
an oversimplified statement. However, the study of single frequency should be helpful for an 
understanding of the very basics of the flow behavior and hence the solute transport in the 
presence of an external magnetic field. Figure 4 contrasts the fluid flow and concentration 
distributions in the melt system with and without the presence of an external magnetic field (Ha 
=120, corresponding to B=0.3 Tesla) at a moment in time after a g-jitter component of g = 1 x 
10' 3 go sin(0.2 n t) sets in. This type of disturbance is known to have a deleterious effect on melt 
flow and hence solutal transport in space. Apparently, even with a moderate magnetic field, the 
melt flow intensity is much reduced, which in turn make$the concentration profile less distorted. 
Analysis of the additional results further illustrates that the period of time harmonic oscillation 
cycle remains the same as the driving g-jitter component, with or without the applied magnetic 
fieldj£ With an applied magnetic field, however, the phase angle is changed and is found to b£ a 
function of the magnitude of the applied field, in addition to other known factors^ such as^ the 
location, the local velocity and the amplitude of the g-jitter component [14]. Also, the 
application of an external magnetic field results in a reduction of transient period over which the 
fluid flow evolves to a quasi-steady state time harmonic oscillation. For pure double diffusive 
convection, it was found that the subharmonic oscillation does not occur until the g-jitter 
perturbations become unrealistically high [14], The applieation of rim, applied magnetic field 
further delays the occurrence of sub-harmonic motions. Sub-harmonic motion comes from the 
nonlinear (or convective) term of the momentum equation. For the conditions being studied, the 
velocity gradient is smaller, although the main flow is strong, thereby causing the flow to 
oscillate with the same frequency as the driving force. Magnetic damping further reduces the 
magnitude of the velocity but the velocity gradient remains approximately^^ same level. As a 
result, subharmonic motion becomes even less likely with an applied field. Numerical 
simulations with different magnetic field strengths also revealed that a stronger magnetic field 
produces a more pronounced damping effect and helps to further reduce the transient period and 
to make the final harmonic oscillation almost in phase with the driving g-jitter force. 

Concentration in the system is strongly affected by the melt convection and the decrease 
of flow intensity should reduce the solutal striation in the melt pool. This is confirmed by Hj^ 
numerical results obtained for the above g-jitter condition with different magnetic fields. Figure 
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5 depicts the transient development of the velocity components and concentrations at a few 
locations as a function of applied fields. Clearly, the use of magnetic field^ produces strong 
damping effect^ on convection and thus smoothes out the solutal fluctuations. Also, a stronger 
magnetic field is more effective in suppressing the local velocity and hence the concentration 
non-uniformity. 


B. Multiple Frequency g-Jitter 

Multi-frequency g-jitter perturbations may be considered as a synthesized representation of 

g-jitter containing a range of frequencies and amplitudes (see Eq. (1)). They are useful for the 

study of the combined effects from various g-jitter components and their interaction with the 

applied magnetic fields. Numerical simulations were made with multiple g-jitter composed of 4 

components with frequencies of (1, 0.1, 0.01, 0.001) and amplitudes of (l.OxlO' 3 , l.OxlQ' 3 , 

1.0x1 O' 4 , l.QxlO''), and with applied magnetic fields. Space experience showed that these g- 

jitter components produce the most detrimental effects on the flow and hence solutal 

distributions in processing systems for space application^" ^Frfen the^Eundamental point of view, 

it is also important to know what effects an applied magnetic field would have on the individual 

\ g-jitter components when they all are present. The results obtained from the numerical model 

,^<06 are' illustrated in Figures 6, 7 and 8. Inspection of these results indicates that the applied 

_ magnetic fields are capable of suppressing the multiple frequency g-jitter induced double 

\-\cj ^ 

diffusive flows and hence are useful in reducing the solutal non-uniformity in the melt pool. 




c ? (pf ( ' Further numerical simulations revealed some additional interesting points that are unique to 

multiple frequency g-jitter with the presence of a magnetic field. First, the applied magnetic 
field is most effective in damping out the influence of the most detrimental g-jitter component, 
which is the component with the largest oscillation amplitude. Early studies [14] on g-jitter 
driven double diffusive convection in the absence of a magnetic field have concluded that for the 
conditions being studied the flow's are essentially a superimposition of the flows by individual g- 
jitter components, and nonlinear effects such as subharmonic motions are a non-occurrence in 
space environment. Also, the component with the largest amplitude generates the most 
intensified convection. This llovv field interacts with the applied magnetic field to give rise\fie 

A 

strongest opposing Lorentz force, thereby producing the most pronounced damping effects. 
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Second, with a large enough g-jitter, which is not realistic but theoretically possible, 
subharmonic motion becomes possible. Such a subharmonic motion can be suppressed by an 
appropriately chosen magnetic field and the nonlinear effect can be completely damped out. 
Moreover, a higher magnetic field produces a more effective damping result, again because of a 
larger Lorentz force. This is the same as-e«e©ttrrtere4 in the case of single frequency g-jitter, as 
discussed above. Finally, solutal non-uniformity decreases with the increasing magnetic field 
strength as a result of higher damping force. With a strong field, the local concentration 
oscillation may be suppressed beyond detection and a diffusion dominant growth condition is 
possible. The precise prediction of such conditions, of course, should have to consider the 
detailed crystal growth configurations and operating parameter, which is beyond the scope of the 
present paper, although the computer code is fully capable of doing so [16]. 


C. Real g-Jitter 

In space environment, crew motions and other space related factors contributed to the 
departure from the perfect dynamic weightlessness condition or g-jitter perturbations. As 
discussed earlier, these perturbations are random in time and in direction (see Figure 1). 
Numerical simulations were performed for the real g-jitter data as given in Figure 1 in 
combination with an applied magnetic field. Some of the results are selectively presented in 
Figures 9, 10 and 11. As one might expect, no time harmonic motion is generated by such 
random g-jitter data. It is seen from these figures that an applied magnetic field is effective in 
reducing the double diffusive convection driven by these random g-jitter forces and hence the 
solutal striation. Several important points may be summarized below for these results. They are 
in principle consistent with the conclusions reached from the single- and multi-frequency studies 
described above. First, the recirculatingg flow pattern is unaltered, but the velocities are reduced 
with an applied magnetic field. Second, the magnetic field is more effective in damping out the 
large component effect or the g-jitter exhibiting the sudden spiking. This is clearly evident in 
Figure 10, where the largest percentage of velocity reduction is associated with the largest g- 
jitter component. Third, the damping effects become more pronounced with the stronger 
magnetic fields. For the conditions studied, the concentration fluctuations caused by spiking in 
g-jitter decrease sharply with the increasing field strength, which is a direct consequence of a 
bigger reduction in velocity. In fact, with FIa=120, the concentration irregularity caused by g- 
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jitter spiking is almost completely smoothed out, as appears in Figure 11. Thus, for the melt 
system, an application of a magnetic field measured by Ha=120 should suppress the melt 
convection ^maH-enough so as to make the solutal transport process a diffusion dominant event. 

5. CONCULDING REMARKS 

This paper has presented a numerical study of the g-jitter driven double diffusive convective 
flow, thermal and. concentration distributions in binary alloy melt systems subject to an external 
magnetic field. The study is based on the finite element solution of the transient 
magnetohydrodynamic equations governing the momentum, thermal and solutal transport in the 
melt pool. Extensive numerical simulations were performed using both the synthesized single- 
and multi-frequency g-jitter as well as the real g-jitter data taken during space flights with or 
without an applied magnetic field. Computed results indicate that the main melt flow is in 
essence a combination of velocity components driven by each individual g-jitter components, 
whether or not an external magnetic field is present. The flow field is characterized by a 
recirculating convection loop, which is caused by both thermal and solutal gradients and 
oscillates in time with a defined frequency equal to that of the driving g-jitter force. An applied 
magnetic field does not alter the recirculating pattern and 'j^or the oscillation period, but it 
changes the phase angle. More important, it interacts with the flow field to produce an opposing 
Lorentz force that in turn suppresses the existing double diffusive convection. As a result, the 
solutal striation and time fluctuations in the melt are reduced. Also, the transient period over 
which the velocity reaches a quasi-steady state time harmonic oscillation is shorter with an 
applied field than with it. The transient period for the concentration field, however, is much 
longer than that for the velocity, and is practically independent of an applied magnetic field. It is 
governed by solutal diffusion, and can be as long as 500 periods of oscillation before the 
concentration field finally reaches a quasi-steady state time harmonic oscillation. For a given 
field strength, the magnetic damping effect is more pronounced on the velocity associated with 
the largest g-jitter component present and/or the g-jitter spiking peaks. Furthermore, damping 
effects (i.e., the reduction of both flow intensity and concentration non-uniformity and 
fluctuations) become more pronounced with increasing magnetic field strength^ A higher 
magnetic field is more helpful in bringing the convection oscillation in phase with the g-jitter 
driving force. Numerical results further indicate that with appropriately selected magnetic field 



strength/ the convective flows caused by g-jitter can be suppressed sufficiently so that the 
solutal transport in the melt is dominated by molecular diffusion. 
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Figure Captions 


Figure 1. Signature of g-jitter data in the x and y directions a function of time registered by an 
accelerometer aboard Space Shuttle, where time is measured in seconds. 

Figure 2. Schematic representation of g-jitter driven oscillating convection in a 2-D cavity in the 
presence of a magnetic field. 

Figure 3. Comparison of finite element solution and analytical solution for a single frequency g- 
jitter induced thermal convection in a parallel plate channel subjected to a transverse magnetic 
field. 

Figure 4. The velocity and concentration field induced by a single frequency g-jitter driving 
force in the y-direction with and without an external magnetic field, (a) Velocity (t=l 4.33 163, 
Ha=0, U max =0.6367), (b) Concentration (t=14.33163, Ha=0), (c) Velocity (t=14.33163, Ha=120, 
U m ax-0.1558) and (d) Concentration (t=14.33163, Ha=120). 

Figure 5. Effects of applied magnetic fields on the time variation of velocity component U, V ' 
and concentration with oscillating gravity, (a, c, e) at (0.5, 0.9220) and (b, d, f) at (0.9220, 0.5). 

Figure 6. The velocity and concentration field under multiple frequency driven force and with 
and without magnetic field, (a) Velocity (t=2. 688364, Ha=0, U max =3.363), (b) Concentration 
(t=2.688364, Ha=0), (c) Velocity (t=2. 688364, Ha=40, U ma x = 0.4488) and (d) Concentration 
(t=2. 688364, Ha=40). 

Figure 7. The dynamic development of velocity and concentration distribution with multiple 
oscillating gravity with and without magnetic field at (0.9220, 0.5): (a, c, e) at Ha=0 and (b, d, f) 
at Ha=120. 

Figure 8. Evolution of concentration as affected by the presence of magnetic field at (0.5, 
0.9220). (a) Ha=0: (b) Ha=40: (c) Ha=120 and (d) Ha=200. 

Figure 9. Magnetic damping effects on the double diffusive convection and solutal distribution 
induced by real g-jitter driven force with and without magnetic field, (a) Velocity (t=0.2698170, 
Ha=0, U max =0.1579), (b) Concentration (t=0.269817(T, Ha=0), (c) Velocity (t=0.269817, 
Ha=200, U, r .ax = 0.2957x 1 0" 1 ) and (d) Concentration (t=0.269817, Ha=200). 

Figure 10. Time variation of velocity component U and V under real g-jitter and with and 
without magnetic field, (a. c) at (0.9220, 0.5) and (b, d) at (0.5, 0.9220). 

Figure 11. Time variation of concentration distribution under real g-jitter and with and without 
magnetic field, (a) concentration at (0.9220, 0.5) and (b) concentration at (0.5, 0.9220). 
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Figure 2. Schematic representation of g-jitter driven oscillating convection with magnetic 
field in a 2-D cavity in the presence of a magnetic field under microgravity condition. 



Figure 4. The velocity and concentration field induced by a single g-jitter driven force in 
y direction with and without an external magnetic field, (a) Velocity (t=T4.33163, Ha=0, 
Umax = 0.6367), (b) Concentration (t=14.33163, Ha=0), (c) Velocity (t=14.33163, Ha=120, 
U m ax-0.1558) and (d) Concentration (t=14. 33163, Ha-120). 
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Figure 6. The velocity and concentration field under multiple frequency driven force and 
with and without magnetic field, (a) Velocity (t=2. 688364, Ha=0, U max = 3.363), (b) 
Concentration (t=2.688364, Ha=0), (c) Velocity (t=2.688364, Ha=40, U max =0.4488) and 
(d) Concentration (t=2. 688364, Ha=40). 






Figure 7. The dynamic development of velocity and concentration distribution with multiple oscillating gravity with and without 
magnetic field at (0.9220, 0.5): (a, c, e) at Ha=0 and (b, d, f) at Ha=120. 
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Figure 9. Magnetic damping effects on the double diffusive convection and solutal 
distribution induced by real g-jitter driven force with and without magnetic field, (a) 
Velocity (t=0.2698170, Ha=0, U ma x=0.1579), (b) Concentration (t=0.2698 1 70, Ha=0), (c) 
Velocity (t=0.269817, Ha=200, U max =0.2957xl0' 4 ) and (d) Concentration (t=0.269817, 
Ha=200). 





